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Abstract 

Lienard systems of the form x + ef(x)x + x = 0, with f(x) an even func- 
tion, are studied in the strongly nonlinear regime (e — > oo). A method 
for obtaining the number, amplitude and loci of the limit cycles of these 
equations is derived. The accuracy of this method is checked in several 
examples. Lins-Melo-Pugh conjecture for the polynomial case is true in 
this regime. 
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1 Introduction 



Many systems in nature display self-sustained oscillations: there is an 
internal balance between amplification and dissipation, and they do not 
require an external periodic forcing to oscillate. For instance, the beating 
of a heart, some chemical reactions, self-excited vibrations in bridges 
and airplane wings, Benard-von Karman vortex street in the wake of 
a cylinder, etc. These phenomena can be modelled by the stable limit 
cycles found in specific nonlinear autonomous dynamical systems. These 
are called 'self-oscillators' 

Limit cycles are isolated closed trajectories in phase space (an inher- 
ently nonlinear phenomenon) 0. They describe the periodic motions of 
the system. A very well known example having one limit cycle is the 
van der Pol equation: x + e(x 2 — l)x + x = 0, where x(t) = dx(t)/dt. 
It displays a wide range of behavior, from weakly nonlinear to strongly 
nonlinear relaxation oscillations when the parameter e is modified, mak- 
ing it a good model for many practical situations || [|, B|. The existence, 
uniqueness and non-algebraicity of its limit cycle has been shown for the 
whole range of the parameter e that controls the nonlinearity || 0, §]. 

Lienard equation, 



with e a real parameter, is a generalization of the van der Pol self- 
oscillator. There are no general results about the existence, number, 




(1) 
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amplitude and loci of the limit cycles of this system || |10|. A nonper- 
turbative method for obtaining information about the number of limit 
cycles and their location in phase space when f(x) is an even polynomial 



is presented in [ 1 1 1 . When / (x) is a polynomial of degree iV = 2n + 1 
or 2n, Lins, Melo and Pugh conjectured (LMP-conjecture) that n is the 
maximum number of limit cycles allowed fl2fl . This conjecture is true if 
f(x) is of degree 2, if f(x) is of degree 3 or if f(x) is even and of degree 
4. Different results about the necessary conditions that certain families 



of f(x) must satisfy to have n limit cycles have been given in |L3| and 
references therein. 

In this work, we are interested in the strongly nonlinear regime (e — > 
oo) of Lienard equation when the viscous term f(x) is a continuos even 
function, otherwise arbitrary. We give a method to find the limit cycles 
in this regime, and we claim that LMP-conjecture is true in that limit. 



2 Lienard Equation 

We start by considering the modified form of Lienard equation (1) after 
the change of variables x(t) = y(x) and x(t) = y(x)y'(x) (where y'(x) = 
dy/dx): 

yy' + ef(x)y + x = 0. (2) 

The variables are now the coordinates (x, x) = (x, y) on the plane. Our 
interest in the shape of limit cycles in phase space (x, y) justifies the 
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elimination of the time variable. 



2.1 Symmetries 

A limit cycle Ci of equation (2) is a closed orbit around the origin (0, 0), 
the only fixed point of the system. This curve C\ = (x, y±(x)) cuts the x- 
axis in two points: (— ai,0) and (02, 0) with 0,1,0,2 > 0. Write (x,y+(x)), 
with y+(x) > and —a\ < x < a 2 , the positive y-branch of the limit 
cycle and (x,y-(x)), with y~(x) < and —a x < x < a 2 , the negative 
one. The inversion symmetry (x,y) <-> (—x, —y) is verified by equation 
(2) because f(x) is an even function. Since the flow lines do not intersect 
themselves, a limit cycle and its transformed by this symmetry must be 
the same curve, then y+(x) = —y^(—x) and a x = a 2 = a. Therefore, in 
the following, we will restrict ourselves to the positive branches of limit 
cycles, (x,y + (x)), with —a < x < a. The amplitude of oscillation will be 
the number a. For a given e this amplitude a identifies the limit cycle 
(Fig. 1). Thus, the number of limit cycles of system (2) is equal to the 
number of different possible amplitudes a. 

Let us remark also the parameter inversion symmetry (e, x, y) <-> 
(— e, x, —y) of equation (2): if C\ = (x,y±(x)) is a limit cycle for a given 
e then Ci = (x,y ± (x)) = (x, —y T (x)) is a limit cycle for — e. Moreover 
if Ci is stable (or unstable) then C\ is unstable (or stable, respectively). 
This is a consequence of the fact that each limit cycle encloses the origin 
(0, 0) and this point changes its stability when e changes of sign. In the 

4 



regime e — > oo the two eigenvalues of the linear part of system (2) in the 
origin are Ai ~ — e _1 /(0) 1 an d ^2 ~ — e/(0). In order to have a well 
defined stability, we impose the condition /(0) 7^ 0. This symmetry is 
used in our computer-simulations to find the unstable limit cycles. For a 
given e the unstable limit cycles are the stable ones for — e after reflection 
on the x-axis. Summarizing, in order to clasify all the limit cycles for a 
given even function f(x), it is enough to find the positive ^/-branch y+(x) 
of the limit cycles of equation (2) when e — > +00. 

2.2 Scaling 

An easier understanding of the behavior of equation (2) in the strongly 
nonlinear regime is obtained by performing the change of variable y = ez. 
With this change Eq. (2) reads: 

zz' + f(x)z = -e- 2 x (3) 

Different scalings can be considered: 
On one hand, if we consider x and f(x) of order 1 in the oscillation region 
—a < x < a of a limit cycle z(x), close to the extreme points (where z is 
of order less than e~ 2 ), the limit cycle obeys the equation: 

zz' + e~ 2 x = yy' + x = 0. (4) 

Integrating this expression we obtain: x 2 + y 2 = a 2 , where a is the 
amplitude of the limit cycle. If x = a + 5x and y = 5y, then close to the 
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extreme points (±a, 0) we have: 2a5x + 5y 2 = 0, with Sy <^ e 2 . 

On the other hand, if z is of order bigger than e~ 2 Eq. (3) is reduced to: 

z(z' + f(x)) = 0. (5) 

It is clear that all relevant information about the limit cycles of Eq. (2) 
when e — > oo is contained in Eq. (5). 



3 Limit Cycle Solutions 



Local solutions of Eq. (5) are: 

Zi(x) = 0, or 
z 2 (x) = -F{x) + C, 

where F(x) = Jq f(t)dt is an odd function and C is a constant. The 
positive |/-branch y + of a limit cycle solution of Eq. (1) and of amplitude 
a will be a solution zi(x) of Eq. (5) verifying (i) z\ (—a) = zi(a) = 
and (ii) zi(0) > 0. The solution Zi(x) = for — a < x < a does not 
represent a limit cycle because it does not verify condition (ii). The 
solution z 2 (x) = —F(x) + C for —a < x < a is not a limit cycle either 
because it does not verify condition (i). Therefore the positive y-branch 
of a limit cycle (solution of Eq. (5)) must be a piecewise function built 
with Zi(x) = and z 2 (x) = —F(x) + C (with C > by condition (ii)) as 
integrant blocks. 
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3.1 Two-Piecewise Solutions 



We can start trying two-piecewise functions called two-piecewise limit 
cycles. Write Xi the gluing point of the limit cycle zi(x) where Z\{x) 
and Z2(x) are glued: Zi(x) and z[{x) must be continous in Xi because the 
velocity x(t) = e _1 ^(x) and the acceleration x(t) = e~ 2 zi(x)z' l (x) are 
continuos in the oscillation. Continuity in zi(x) means C = F(xi) > 0. 
Continuity in z[{x) means /(xj) = 0. The point — is also a zero of f(x) 
(in the following we will suppose f(x) has a finite number of zeroes) and 
it is the gluing point of the negative y-branch y_ of the limit cycle. Thus 
each pair of zeroes ±Xj of f(x) can generate at most one limit cycle. If 
f(x) is a polynomial of degree 2n, there will be at most n (two-piecewise) 
limit cycles. 

The stability of a limit cycle y(x) solution of Eq. (2), with —a < x < 
a, is determined by the sign of the integral: 

J-a y(x) J-a Z(X) 

If a < the limit cycle is stable and if a > it is unstable. In the 
two-piecewise case this quantity a is controlled essentially by the region 
where z(x) = 0, more exactly by the sign of f(x) close to the gluing point 

= -sign[f(x ~ Xi)} z(x)=0 . 

z(x)=0 

In the next section it is shown that if Xi is a gluing point, F(xi) must be 
a maximum of F(x). Therefore, for e > 0, if X{ = Sj < the limit cycle 
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sign [a\ 



-sign 



z(x) 



dx 



whose positive y-branch is y+(x) = ezi(x) is stable because, at the side 
where z(x) = 0, f(x ~ Sj) > 0. If Xi = Ui > it is unstable because, at 
the side where z(x) = 0, /(x ~ Uj) < 0. 



Thus the expression of the positive y-branch = €zf(x) (solution 

of Eq. (5)), with — af < x < af, of a stable (two-piecewise) limit cycle is 
(Fig. 2(a)): 

s( \ _ J if - af < x < Si , . 

W ~ \ -F(x) + F( Si ) if < rr < W 

and the form for an unstable (two-piecewise) limit cycle, y+(x) = ezf(x), 
with -a" < x < is (Fig. 2(b)): 

, , _ f -F(x) + F( Ui ) tf-a?<x< Ui , 
z ^ x )-\ if Ul <x<a^ K ' } 

where af' u > represent the amplitude of each limit cycle, and and U; t 
are the gluing points of the two pieces of each cycle, zf and zf, respec- 
tively. Recall that Sj < for the stable cycle and > for the unstable 
one. In both cases F(si) > 0, F(ui) > and /(si) = /(««) = 0. 



3.2 Number of Two-Piecewise Solutions 



STABLE CYCLES : We study in more detail the stable limit cycles given 
by Eq. (6). The coordinate y+(x) = ezf(x) of the limit cycle vanishes in 
the extreme points ±af and then F(si) = F(af). As zf(x) is positive for 
Si < x < af the amplitude af is defined by: 

af = min {x > Sj, F(x) = F(si)} . 
8 



Also the property af > \s{\ is fulfilled. 

Let us invert the problem. If s* < verifies f(s*) = and F(s*) > 0, 
is it possible to build a stable (two-piecewise) limit cycle, as given by Eq. 
(6), with s* as gluing point?. First we define a* by the rule: 

a* = min {x > s*, F(x) = F(s*)} . 

Geometrically a* represents the x-coordinate of the first crossing point 
between the straight z = F(s*) and the curve F(x) (Fig. 3(a)). If 
a* < \s*\ it is not possible to build the limit cycle and we can eliminate 
this s* as a possible gluing point. If a* > \s*\ the point s* is a gluing point 
candidate. We rename all the pairs (s*,a*) verifying this last property 
as (si,aj) and collect them into the set: 

A s = {fa, a?), f(si) = 0, Ffa) > 0, s i+1 < s t < 0, a* > \si\} . (8) 

By construction af +1 > a\. 

There are two different situations when two sucessive pairs, (s«, af) 
and (s i+ i,a s i+1 ), are ordered: 

(a) — a s i+1 < Sj+i < —af < Sj. In this case it is possible to build a two- 
piecewise limit cycle with the pair (si,af) as indicated by Eq. (6). This 
pair is picked out and renamed once more as (s^ af). 

(b) — a| +1 < — af < Sj+i < Sj. Now the constuction of a limit cycle 
derived from the pair (sj,af) is not possible. If a initial condition (with 
z = 0) in the interval [—af, s i+1 ] is given, the system will jump from the 
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point Sj+i to af +1 > af, then the curve does not close at af and there is 
no a limit cycle with amplitude af. This pair is rejected. 

If there is only one pair (s±, a{), we consider it satisfies (a). 

All the existing (two-piecewise) stable limit cycles can be found com- 
paring the pairs i and i + 1 under rules (a)-(b) and iterating this process. 
All the pairs selected by condition (a) (and renamed as (s i: af)) are col- 
lected into the set: 

A s = {( Si , af)} = {{s h af) G A% (s h af) verifies (a)} . (9) 

The algorithm above proposed shows also that F(x) must be a local 
maximum at each gluing point Sj G A s . 

The number, l s = card(A s ), of pairs (sj, af) is the number of stable 
(two-piecewise) limit cycles of system (5). 

UNSTABLE CYCLES : The same process can be repeated for the 
unstable cycles by considering the points u* > 0, where f(u*) = and 
F(u*) > and finding their partners a* defined as (Fig. 3(b)): 

a* = max {x < u*, F(x) = F(u*)} . 

Moreover the gluing point candidates must verify \a*\ > u*. After col- 
lecting the pairs (u*,a*) fulfilling this last condition, we have, as in (8), 
the set: 

A u = {(Hi, <), f(Ui) = 0, Ffa) > 0, u i+1 >Ui>0, K| > Ui] . (10) 
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A similar algorithm as indicated above can be applied in this case with 
the following modified rules: 

(a') ui < —af < u i+ i < —af +1 . In this case there exists an unstable 
two-piecewise limit cycle resulting from the pair (tZj, af) and given by 
Eq. (7). This pair is picked out and renamed (ui,af). 
(b') U{ < u i+ i < —af < —af +1 . The pair (-Uj, af) does not produce a such 
limit cycle and is rejected. 

If there is only one pair (ui,af) we consider it satisfies (a'). 

We iterate the process given by rules (a')-(b'). All the pairs selected 
by condition (a') are collected into the set: 

A u = {(u h af)} = {(u h af) G A u , (u h af) verifies (a')} . (11) 

The algorithm above proposed shows also that F(x) must be a local 
maximum at each gluing point G A u . 

The number, /„ = card(A u ), of pairs (ui,af) is the number of unstable 
(two-piecewise) limit cycles of system (5). Obviously, l s — \<l u <l s + \. 

We claim that the total number / of limit cycles of Eq. (1) in the 
strongly nonlinear regime is I = l s + l u , where l s and l u are the number of 
stable and unstable limit cycles of Eq. (1), respectively. The amplitudes 
of these limit cycles are given by the numbers af and af, respectively. 
We remark also that each pair of zeroes ±Xi of f(x) produces at most one 
limit cycle. If f(x) is a polynomial of degree 2n there will be at most n 
limit cycles. Therefore, LMP-conjecture is true in the strongly nonlinear 
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regime. 



3.3 Shape of the Limit Cycles 

The method introduced above allows us to find the number / of limit 
cycles and their amplitudes. However, two-piecewise solutions (Eqs. (6)- 
(7)) are only a first approach to the problem of finding the exact loci of the 
limit cycles of equation (1) in the e — > oo regime. A better approximation 
to the shape of the limit cycles is presented in this section. 

STABLE CYCLES : A stable two-piecewise solution zf(x) of Eq. (5), 
identified by the pair (sj, af), is built of two blocks glued at sf z^ { (x) = 
and Z2i(x) = —F(x) + F(si) (see Eq.(6)). The role of the point s$ is to 
allow the limit cycle to jump from the left side where zf(x) = zf^x) = 
to the maximal amplitude af through zf(x) = z 2t i(x). Depending of f(x), 
it is also possible to find some points inside the interval — af < x < 
and with similar jumping properties than Sj. If this is the case, the exact 
limit cycle in the interval [— af,Sj] is not zf(x) = z{^(x) = 0, and the 
piece zl ^x) in this interval must be corrected in the following way (Fig. 
4): 

(a) Define {s^i, • • • , Sj jfli }, the largest set of points where F(x) is a local 
maximum and verifying: 

f(si,i) = ■■■ = f{si, ni ) = 0, 
-a- < Si,i < s i)2 < • • • < s ij7lj < ^, 
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< F(s itl ) < F(s it2 ) <■■■< F(s itni ) < F( Si ). 



The values • • • , Sj inj } will be called the jumping points of the cycle 
(si,af). 

(b) For each pair (i,m), m = l,---,nj, we define bi )Tn = min{x > 
Si, m ,F(x) = F(si tm )}. It follows that bi tm < Si jTn+ i for every m because 
F(s ijm ) < F(si >m+ i). Thus we have the ordered set: 

-Oj < Si,i < < s ij2 < 6i,2 < • • • < s i)Ui < b ijni < Si. 

(c) The correct expression zf(x) of the stable limit cycle (af , Sj) solution 
of Eq. (5) is: 



where 





if 


— af < x < s,i 


~ | —F(x) + F(sj) 


if 


Si < x < af, 





if 


— af < x < Si t i 




if 


S^i < x < b it i 





if 


b it i < X < Sj, 2 


-F(x) + F(s i>2 ) 


if 


Sj, 2 < X < foj )2 





if 


6 i>2 < a; < Sj i3 


-F(x) + F( Si , ni ) 


if 


^i,n, *^ ^ ^i,r 





if 


^i,nj < 2? < Sj. 



(12) 



Note that the dynamics in the jumping points s^ m has the correct 
scale to be governed by Eq. (5). Then the continuity in z(x) and z'(x) 
can be imposed in these points in the gluing process. Nevertheless, a 
more subtle study precises this gluing process in the points b^ m because 
other scales not present in Eq. (5) take place in the dynamics. Thus 
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the continuity of z'{x) (the acceleration) is lost in the points b^ m . This 
fact is not a problem in order to consider the expression (12) as a fine 
approximation (up to order e~ 2 ) to the shape of the stable limit cycle, 
y+(x) = ezf(x), solution of Eq. (1) in the strongly nonlinear regime. 

UNSTABLE CYCLES : Similarly to the previous case, the piece z^x) 
in the unstable two-piecewise solution (7) must be corrected. If the limit 
cycle is identified by the pair (ui, af) the correction in the interval [v,i, a"] 
is given by the following steps (Fig. 4): 

(a) Define {u it i, • • • , u itPi }, the largest set of points where F(x) is a local 
maximum and verifying: 

f(ui,i) = ■■■ = f(u i<Pi ) = 0, 
Ui < u i:1 < u i:2 < ••• < u iyPi < a", 
F{ Ui ) > F(ui tl ) > F(u h2 ) >■■■> F(u itPi ) > 0. 

The values {u^i, • • • , Ui tPi } are the jumping points of the cycle (u^ a"). 

(b) For each pair (i,m), m = 1,- ••,£>,, we define d i>m = max{i < 
u i>m ,F(x) = F(u itm )}. It follows that d^ m > u itm _i because F(u i>m ) < 
F(ui tm -i). Thus we have the ordered set: 

Ui < d it i < u it i < d i)2 < u i)2 < ■■■ < d i)Pi < u ijPi < a". 

(c) The correct expression zf{x) of the unstable limit cycle (a", Wj) solu- 
tion of Eq. (5) is: 

_ f -F{x) + F( Ui ) tf-a?<x<Ui 
z ^ x >-\ z v.( x ) if Ui<x <a?, [6) 
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where 



if 

F{x)+F{u hl ) if 

if 

F{x)+F{u h2 ) if 

if 



Ui < x < d it i 
d it i < x < u iA 
Ui t i < x < di i2 
d ij2 < x < u i:2 
u i>2 < x < d ij3 



F(x)+F(u ijH ) if 
if 




Similar comments to those of the stable case apply here. Therefore 
the expression (13) is a fine approximation (up to order e -2 ) to the shape 
of the unstable limit cycle, y+(x) = ezf(x), solution of Eq. (1) in the 
strongly nonlinear regime. 

4 Examples 

We illustrate in this section the method introduced in Sections 2-3 for 
finding the number, amplitude and loci of the limit cycles of equation (1) 
by means of some examples. 

(1) f(x) = 5x 4 — 3x 2 — 1 , then F(x) = x 5 - x 3 - x. The only real 
solutions of f(x) = are x = ±0.9157. The only local maximum verifying 
F(x) > is s* = -0.9157. The only value a* verifying F(a*) = F(s*) is 
a* = 1.3837. Moreover, it verifies a* > |s*|, then we can rename this pair 
as (s*, a*) = (si, di) and because there is only one pair, (s 1 , d{) = (si, a\). 
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The limit cycle when e — > oo is then (F(si) = 1.0397): 
zi(x) = 



if - 1.3857 < x < -0.9157 

x 5 + x 3 + x + 1.0397 if - 0.9157 < x < 1.3857 



The amplitudes a exp for different values of the parameter e can be calcu- 
lated by numerical integration of Eq. (1). For instance: 



e 


^exp 


1 

10 
100 


1.4099 
1.3975 
1.3874 



Let us remark the agreement between our analytical approach and the 
behavior of the system. 



(2) f(x) = 5x 4 - 3x 2 + 0.1 , then F(x) = :r 5 -:r 3 ±0.1:r. The solutions 
of f(x) = are x 1 = ±0.7514 and x 2 = ±0.1882. The only local maxima 
verifying F(x) > are s* = -0.9157 and u* = 0.1882. In the case of s* 
the only value a* verifying F(a*) = F(s*) is a* = 1.0045. Moreover, it 
verifies a* > \s*\, then we can rename this pair as (s*,a*) = (si,al) = 
(si,a{). The stable limit cycle, when e — > oo, is then (-F(si) = 0.1096): 

/ \ _ ( if - 1.0045 < x < -0.7514 

Zl( - X) ~ \ -x 5 ±x 3 - 0.1x + 0.1096 if - 0.7514 < x < 1.0045 

The amplitudes a exp for some values of the parameter e are: 



e 


&exp 


1 

10 
100 


1.0234 
1.0164 
1.0096 
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In the case of u* the closest point a* < to the origin verifying F(a*) = 
F(u*) is a* = —0.3945. Moreover, it verifies \a*\ > u*, then we can 
rename this pair as (u*,a*) = (ui,af) = (si,a"). The unstable limit 
cycle when e — > oo is then (F(ui) = 0.0124): 

f -x 5 + x 3 - O.lx + 0.0124 if - 0.3945 < x < 0.1882 
Zl[ - X > ~\ if 0.1882 < x < 0.3945 

The experimental amplitudes a exp for different values of the parameter e 
are: 



e 


^exp 


1 

10 
100 


0.3909 
0.3927 
0.3943 



(3) f{x) = -{x 2 - 0.09){x 2 - 0A9){x 2 - 0.81) , thenF(a;) = -0.1428x 7 + 
0.278x 5 - 0.1713a; 3 + 0.3572a;. The solutions of the equation f(x) = 
are: x = ±0.3, x = ±0.7 and x = ±0.9. The only local maxima are 
«* = 0.9 and u* p = 0.3. The value a* = 0.5443 verifies F(a* a ) = F(u* a ) 
and | a* | < m* = 0.9 then this pair can not build a limit cycle. On the 
other hand, the value a*p = 1.0485 is the closest point to the origin verify- 
ing F(a*p) = F(u%). Moreover, it verifies \a%\ > u%. Then, we can rename 
this pair (u%,a%) = (-Ui,a") = (ui,a^). As u\ < |a"| the value u\ is the 
jumping point u± t i = 0.9 of the cycle {u\,a\). The only (unstable) limit 
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cycle of this system is therefore (F(ui) = 0.0067 and F(ui t i) = 0.0031): 



zi{x) 



0.1428x 7 - 0.278x 5 + 

0.1713a; 3 - 0.3572a; + 0.0067 if 

if 
0.1428a; 7 - 0.278a; 5 + 

0.1713a; 3 - 0.3572a; + 0.0031 if 

if 



- 1.0485 < x < 0.3 
0.3 < x < 0.5433 

0.5433 < x < 0.9 
0.9 < x < 1.0485 



The experimental amplitudes a exp obtained by direct integration are: 



e 




500 

1000 

5000 


1.0497 
1.0493 
1.0489 



5 Conclusions 



Periodic self-oscillations can arise in nonlinear systems. These are repre- 
sented by isolated closed curves in phase space that we call limit cycles. 
The knowledge of the number, amplitude and loci of these solutions in a 
general nonlinear system is an unsolved problem. 

In this work, we have studied the Lienard equation in the strongly 
nonlinear regime. An effective algorithm for obtaining its limit cycles 
solutions (number, amplitude and loci) has been proposed. There exists 
an strong agreement between our analytical approach and the numerical 
integration of the system. Moreover, we claim that Lins-Melo-Pugh con- 
jecture is true in this regime when the nonlinear viscous term f(x) is an 
even function. 
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Figure Captions 



Figure 1: A typical limit cycle C\ = (x,y±(x)) of Eq. (1). The 
oscillatory dynamics x(t) given by this solution verifies —a < x < a, 
where a is the amplitude of oscillation. Let us observe the inversion 
symmetry (x,y) <-> (—x, —y) of this curve. 

Figure 2: Diagrams (a)-(b) represent two-piecewise limit cycle so- 
lutions of Eq. (5), a stable and an unstable, respectively. Its amplitudes 
are af and a". The functions z\{x) and z 2 (x) = —F(x) + C are glued at 
the points ±Sj in the stable case (a) and at the point ±Uj in the unstable 
case (b). 

Figure 3: Diagrams (a)-(b) show algorithms for obtaining the limit 
cycles, stable and unstable, respectively, (a) The pair (s*,a*) is rejected 
because a* < \s*\. A s = {(si,a?),i = 1,2,3}. A s = {{s 3 ,a s 3 )}. (b) The 
pair (u*,a*) is rejected because \a*\ < u*. A u = {(ui,af),i — 1,2,3}. 
A u = {(ui,af),i = 2,3}. 

Figure 4: Shape of the limit cycles, (a) There are four limit cycles: 
two stables collected into A s = {(s i: af) : i = 1,2,}, and two unstables 
into A u = {(ui,af),i — 1,2}. (b) The unstable cycle {ui,af) has the 
jumping point u 1} \. The stable one (52,^2) nas ^ ne jumping point s 2 ,i. 
Let us remark the repeating shape of the four limit cycles. 
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